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Abstract. Here we prove that Benford's law holds for coefficients of an infinite class 
of modular forms. Expanding the work of Bringmann and Ono on exact formulas 
for harmonic Maass forms, we derive the necessary asymptotics. This implies that the 
unrestricted partition function p(n), as well as other natural partition functions, satisfy 
Benford's law. 



1. Introduction and Statement of Results 

It has long been observed that many naturally occurring statistics and arithmetic 
functions have some surprising properties. In 1881, astronomer Simon Newcomb no- 
ticed that the earlier pages in logarithm tables were more worn than later ones. Instead 
of the a priori estimate that the first digit will be 1 one-ninth of the time, he found in 
these instances that this frequency is approximately 30% for the digit 1 and less than 
5% for the digit 9 [TU]. This phenomenon, known as Benford's Law, appears in a wide 
class of data including river lengths and population demographics. For a more detailed 
discussion on the history and previous work on the subject see PQ, [6], [7]. 

Although this "law" is well-known, it has only been proven to hold for a relatively 
small class of arithmetic functions. For example, Miller and Kontorovich prove Benford's 
law for distributions of values of //-functions and the 3x + 1 problem in [8] . The purpose 
of this paper is to prove that the statistically observed frequencies dictated by Benford's 
law hold for an infinite class of sequences arising as the coefficients of modular forms, 
including the partition function p(n). 

We define a partition of a non-negative integer n 6 N to be any non-increasing se- 
quence of positive integers which sum to n. The partition number p(n) is the number of 
partitions of n. For a sequence of positive integers a(n), let 

. , ... # { n < x : first digits of a(n) in base k are the string d\ 
(1) B(d,x,k;a(n)) = . 
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We say a{n) is Benford if lim B(d,x,k;a(n)) = (\og k (d + 1) — \og k (d)) (mod 1) for 

all integers k > 2. We denote the space of Benford functions as B. Note that for a 
function to belong to B, Benford's Law must hold for any initial string of digits in any 
base. We include this level of generality to study frequencies such as 5(101, x, 2; a{n)) 
which counts the proportion of n < x for which a(n) begins with the string <i="101" in 
base 2. Since 101 2 = 5i and log 2 (6) — log 2 (5) ~ 0.263 (mod 1), the predicted frequency 
for members of a Benford sequence to begin with 101 2 is about 26.3%. 

The following data for initial digits illustrates the plausibility of Benford's Law for 
p(n) for k = 10. 

Table 1. B(d,x, 10;p(n)) 



x d = 1 2 3 4 5 6 7 8 9 

W 033 0T6 Oil O09 O07 0l)6 O07 O05 0.03 

10 3 0.305 0.177 0.127 0.094 0.076 0.068 0.057 0.052 0.044 

10 4 0.302 0.177 0.126 0.096 0.078 0.067 0.057 0.051 0.046 

oo? 0.301 0.176 0.125 0.097 0.079 0.067 0.057 0.051 0.046 



We also have the following data for the first three digits of p(n) in base 2. 





Table 2. 


. B(d, 


x,2;p(n)) 


X 


d = 100 


d=101 


d=110 


d=lll 


200 


0.285 


0.270 


0.205 


0.225 


400 


0.308 


0.273 


0.209 


0.205 


600 


0.313 


0.267 


0.217 


0.198 


800 


0.314 


0.263 


0.219 


0.201 


1000 


0.315 


0.262 


0.220 


0.200 


5000 


0.321 


0.264 


0.222 


0.194 


; 


I 


; 


4- 


; 


oo? 


0.322 


0.263 


0.222 


0.192 



To start, we begin with the following definition and theorem. 
Definition. We say that an integer-valued function a{n) is good whenever 

a(n) ~ b{n)e c ^ 

(where f(x) ~ g(x) means that lim f(x)/g(x) = 1) and the following conditions are 

x— >oo 

satisfied: 

(1) There exists some integer h > 1 such that c(n) is /i-differentiable and c^ h \n) 
tends to zero monotonically for sufficiently large n. 

(2) lim n\c (h \n)\ = oo 

n— >oo 

]j{ h ) logfe(n) 

(3) lim ,, . . . = 0, where denotes the h th derivative. 

V ' n^oo c O)(n) 
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Our first result is the following theorem. 
Theorem 1.1. If a(n) is good, then a(n) G B. 

As a special case of this result we obtain the following corollary. 
Corollary 1.2. The partition function p(n) G B. 
Proof. Using the celebrated Hardy-Ramanujan asymptotic 

it immediately follows that p(n) is good and hence Benford. □ 

Next, we explicitly demonstrate that a large class of arithmetic functions arising from 
the coefficients of modular forms is Benford, as in the following theorem. 

Theorem 1.3. Suppose that M(z) is a weakly holomorphic modular form (see §3.1 for 
the definition) of weight ^ > k G with integral Fourier coefficients and at least one 
pole. Then the nonzero coefficients of M(z) are Benford. 

Generalizing Corollary II. 2\ we obtain the following corollary. 

Corollary 1.4. For any positive integer we say that a partition is s-regular if it has no 
part divisible by s. Denote by b s (n) the number of s-regular partitions. Then Theorem 
11.31 implies that b s (n) is Benford. 

Applying Theorem 11.11 to generalized Dedekind Eta-products, we have the following 
corollary. 

Corollary 1.5. For 5 > 2 and < g < L^jrJ define r g ^(n) to be the number of partitions 
of n into parts congruent to ±g (mod 5). This includes the famous Roger s-Ramanuj an 
functions 



oo 



qTi +an — 



n-l 



for a = 0, 1, where for n > 0, (a; q) n := J^[(l — a><f) is the g-Pochhammer symbol. Note 



i=0 



that these equalities are given by the celebrated Roger s-Ramanuj an identities. Then 
r gt s(n) is Benford. 

Remark. 

It is well-known that every nonconstant weakly holomorphic modular form with non- 
positive weight has a pole. Though this is not true for weights k > \, Theorem 11.31 
applies to all weakly holomorphic modular forms of weight | with a pole and essentially 
all weakly holomorphic modular forms with higher weight which have a pole. See the 
discussion and example in §3.2 for more details. For a more detailed discussion on weakly 
holomorphic modular forms see §3.1 and jllj . 
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2. Theorems on Uniform Distribution 

2.1. Preliminaries on Uniform Distribution and Properties of Benford Spaces. 

We first note that proving a sequence is Benford reduces to a problem of uniform dis- 
tribution. Using this formulation and some classical theorems on uniform distribuiton, 
we derive a set of sufficient conditions for being Benford. 

Definition. For a real sequence a(n), JVeN, and £CR, let A(E, N,a(n)) := < 
iV | a(n) G E}. Then a(n) is said to be uniformly distributed mod 1, if, for all intervals 
[a, b) C [0,1), 

A([a,b),n,{a(n)}) 
hm = o — a, 

where < x > denotes the fractional part of x. 

We begin by recalling a result of Diaconis that a(n) G 23 if and only if log fc (a(n)) is 
uniformly distributed mod 1 for all [4J. To prove our results we need the following 
preliminary theorem. 

Theorem 2.1 (Weyl's Criterion). The sequence a(n) is uniformly distributed mod 1 if 
and only if 

N 

v— x Y 

n=l 

for any / fc 6 Z. 



1 JV 

Jim ±y e 2^.a(n) =Q 



As a corollary to Van der Corput's Difference Theorem we recall the following theorem: 

Theorem 2.2 (Theorem 3.5 of [9|). Let ieN and f(x) be a function defined for x > 1 
which is k-times differentiable for all x > x$ for some xq E M+. Suppose that f^ is 
eventually monotonic, 

lim f {k) (x) = 0, 

x— >oo 

and 

lim x|/ (fc) (x)| = oo. 

x— >oo 

T/ien {f{n) : n G N} vis uniformly distributed mod 1. 
Finally, we note a few basic facts regarding Benford spaces. 



Lemma. 

(1) // f(n) eB andceR- {0}, then cf(n) E B. 

(2) If f{n) G B is nonzero, then G B. 

(3) If f{n) G B and f(n) ~ #(n) ; i/ien gr(n) G 23. 

Note that (1) and (2) follow directly from the definition of uniform distribution and 
(3) follows easily using Weyl's Criterion. 
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2.2. Proof of Theorem ll.il , Here we prove Theorem 11.11 using these results on uniform 
distribution. 

Proof. First note that by Part 3 of the lemma, it suffices to show that b(n)e c ( n > is Benford. 
This is equivalent to showing the uniform distribution of log(6(n)) + c(n) (mod 1). We 
see that c(n) satisfies the limit conditions of Theorem 12.21 by definition. By the order 
assumption on D^ h ' (log(fo(n))) , the limits are unaffected upon adding logfe(n). □ 

3. Proof of Theorem 11.31 

3.1. Weakly Holomorphic Modular Forms. Here we give some of the preliminary 
definitions and theorems on weakly holomorphic modular forms. For a more detailed 
reference, see [TTJ. Throughout this discussion let q := e 2mz . 

For any positive integer N, we define the level N congruence subgroups o/SL 2 (Z) by 

r (iV) := { G SL 2 (Z) :cEOmodN} 



r x (N) := { ) G SL 2 (Z) : a = d = 1 mod N and c = mod n} 

T(N) := | Q ^ G SL 2 (Z) : a = d = 1 mod and 6 = c = mod n} 

Suppose that T is a congruence subgroup of SL 2 (Z). A cusp is an equivalence class of 
P 1 (0) under the usual fractional linear action of T. 



Let f(z) be a meromorphic function on the upper half plane EI, k G Z, and T be a 
congruence subgroup of level N. First, we define the "slash" operator by 

U\ kl ){z) :=(det 7 ) fc / 2 (^ + d)- fc /(7^), 

where 7 G SL 2 (Z) and 72 := jg* 

Then /(z) is said to be a meromorphic modular form with integer weight k on T if the 
following hold: 

'a 6 N 



(1) For all 2 G HI and all ^ ^ G T we have that 

(2) If 70 G SL 2 (Z), then (/|fc7o)(^) has a Fourier expansion of the form 

(/l*7o)W= E a 7oH^- 



If = 0, then /(z) is known as a modular function on T. 

We say that /(z) is a weakly holomorphic modular form if its poles are supported at 
the cusps of T. 

Here we recall the notion of a modular form of half- integral weight. First we define 
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and e d . If d is an odd prime, let f|J be the usual Legendre symbol. For positive 
odd d, define f |J by multiplicativity. For negative odd d, we let 




if d < and c> 0, 



d J ifrf<0andc<0. 



Also let (grr) = 1- Define e^, for odd d, by 



1 if d = 1 mod 4, 
z if d = 3 mod 4. 



Let A be a nonnegative integer and A" a positive integer. Furthermore, suppose that 
X is a Dirichlet character modulo 4A^. A meromorphic function g(z) on EI is said to be 
a modular form with Nebentypus x an d weight X + \ if it is meromorphic at the cusps 
of r, and if 

for all (« *) 6 r (4W). 

Remark. As in the integral case, we say that a meromorphic modular form M(z) is 
weakly holomorphic if its poles are supported on the cusps (that is, M is holomorphic 
on H). 

For example, consider Dedekind's eta-function, 

oo 

rj(z) :=q 1/2i l[(l-q n ), 

n=l 

which is a non- vanishing, weakly holomorphic modular form with weight 1/2. The 
following theorem becomes useful in generating more examples of modular forms (see 
for example, [TTj). 

Theorem 3.1. Let f(z) = YIs\n v(^ z Y s be an eta-quotient with k = §$2<5|iv r <5 G Z. 
Then if 

^25r 5 = (mod 24) 

5\N 



and 



then f{z) satisfies 



EN 
—r 8 = (mod 24), 
o 

S\N 



<^±h j =x{d){cz + d) k f{z) 
cz + dJ 



BENFORD'S LAW FOR COEFFICIENTS OF MODULAR FORMS AND PARTITION FUNCTIONS 7 



for every ^ G T (N). Here the character x is defined by x(d) := l \ s ^ ), where 

3.2. Proof of Theorem I1.3L We are now in position to prove Theorem II .31 

Proof. We first consider the case where k is nonpositive. By hypothesis, M has at least 
one pole at a cusp of the relevant subgroup T C SL^Z). Recall that a harmonic Maass 
form is a function defined in a similar manner as a modular form, but with relaxed 
growth conditions and the additional requirement that it lies in the kernel of a certain 
weight k Laplacian operator. Every such harmonic Maass form can be decomposed as 
a sum of a holomorphic part f^ and a nonholomorphic part It is well-known 

that every weakly holomorphic modular form is a harmonic Maass form. In addition, the 
nonholomorphic part f(~' is identically zero. By a recent paper of Bringmann and Ono 
[3] , we have exact formulas for the holomorphic part of harmonic Maass forms for k < | . 
Because in our setting f(~' is identically zero, these give us exact formulas for the Fourier 
coefficients of M. The explicit formulas are quite complicated, so for brevity note that 
if we write M := Y2 n >n m { n ) ( l n i ft is easy to show that there exists an integer t > 1 
such that for all < r < t the nonzero coefficients satisfy 



m(tn + r) ~ K(M, r, t) ■ (tn + r)~ ■ 7i_ fc (a<(M, r, t)y/tn + r) 

where K(M, r, t) and a(M, r, t) are constants and I denotes the modified Bessel function 
of the first kind. Using the standard asymptotics for Bessel functions, namely that 

we see that the sequence of nonzero coefficients is good and hence Benford. 

If k = 1/2, then by Theorem 1.2 in [3J, we have a new exact formula of the same form 
with the addition of a finite number of terms which do not affect the asymptotic. □ 

Remark. Although we have stated and proved Theorem 1.3 for weights k < |, it turns 
out that a suitably modified version holds for all weights. Let M(z) be a weakly holo- 
morphic modular form with weight k > | with integral coefficients and at least one pole, 
which by definition must be at a cusp. We can decompose M(z) = P(z) + H(z), where 
P{z) is a linear combination of Maass-Poincare series and H (z) is a holomorphic mod- 
ular form. Note that by [3] this decomposition is stated only for the case 1/2 > k e |Z, 
but in fact it will hold for all k e |Z. When k < 0, it follows that H(z) = 0, and 
when k = 1/2, H(z) is a linear combination of readily understood theta functions by 
the Serre-Stark basis theorem [3] . Therefore, when k < | this decomposition gives exact 
formulas for the coefficients. For other k, we do not have a good theory for the possible 
H(z), therefore we typically cannot obtain exact formulas for the coefficients of M(z). 
However, it is well-known that the coefficients of all such H(z) are bounded by a fixed 
power of n. To see this, note that the space of holomorphic modular forms is spanned 
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by Eisenstein series and cusp forms. From the expansion of Eisenstein series in terms 
of divisor functions, it is easily seen that the coefficients are bounded by a polynomial. 
Cusp forms are also bounded by a power of n by Proposition 8 on page 23 of [13] . Con- 
sequently, we obtain sufficient asymptotics for all coefficients which are nonvanishing 
in P(z). These asymptotics will have subexponential growth, which will dominate the 
coefficients from H(z). 

Example. To illustrate the above remark, consider the weight 12 weakly holomorphic 
modular form M(z) := j(z) ■ E^z) 3 , where j(z) is the j-modular invariant and E^{z) is 
the Eisenstein series of order 4. The first few terms are 

M(z) = q' 1 + 1464 + 911844g + 313589120g 2 H 

This function has a pole at oo so the above remark applies as is illustrated in the base 
3 case below. 



Table 3. B(d,x,3; j(z) ■ E A (zf) 



X 


d = 10 


d=ll 


d=12 


d=20 


d=21 


d=22 


500 


.2440 


.2020 


.1700 


.1320 


.1300 


.1220 


1000 


.2590 


.2050 


.1610 


.1360 


.1270 


.1120 


1500 


.2627 


.2027 


.1633 


.1373 


.1273 


.1067 


2000 


.2635 


.2030 


.1660 


.1375 


.1245 


.1055 


1 


I 


1 


I 


I 


1 


1 


oo? 


.2619 


.2031 


.1660 


.1403 


.1215 


.1072 



4. Proof of Corollaries 1.4 and 1.5 

In this section we use Theorem 1.1 to prove new examples of Benford sequences in 
interesting cases. 

4.1. Proof of Corollary 1.4. First note that 

n=0 n=l H 

Then in terms of Dedekind's eta function, i](z), we have that 

,24n+ S -l _ V(2ASZ) 



^b s {n)q 2 



n=0 



7](2Az) 



Then by Theorem 1.64 [11] , we see that b s (n) is obtained as the coeffients of a non- 
constant modular form of weight zero on r (576s). Moreover, as r](z) is nonvanishing 
on H, the given modular form is weakly holomorphic and hence b s (n) is Benford. 
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4.2. Proof of Corollary 1.5. It is easy to see that 

oo 1 

5>*«w= n j]—^) n a 



n=0 l<n=g (mod 5) v y l<«=-3 (mod <5) v y 

For example, we have the Roger s-Ramanuj an functions 

gn 2 +an °° 



„_o fa ?)» n A i (i - g 5n+a+1 )(i - q 5n+4 - a ) 

for a = 0, 1. 

Our claim is that for all g and 5 satisfying the above restrictions r g> s is Benford. To 
see this, define the generalized Dedekind Eta-product T)g,s(z) by 

Vg ,s(z) := e 2 ^ ^ J] (l-q n ) J] (1 - q n ). 

l<n=g (mod S) l<n=—g (mod <5) 

Here P2{t) is the second Bernoulli function P2(t) := {t} 2 — {t} + I, and {t} denotes 
the fractional part of t. Then by [12], for 7^ 0, |5, it follows that T]g,s( z ) is a modular 
form of weight zero which is nonvanishing on H. The generating functions for r 9j< 5 are 
then given by the Fourier expansion of ■ — shifted by an integral power of q. As this 

is a nonconstant weakly holomorphic modular form of weight zero, it follows that r g ^(n) 
is Benford. 
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